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ABSTRACT 

Aims. For many years feedback processes generated by OB-stars in molecular clouds, including expanding ionization fronts, stellar winds, or 
UV-radiation, have been proposed to trigger subsequent star formation. However, hydrodynamic models including radiation and gravity show that 
UV-illumination has little or no impact on the global dynamical evolution of the cloud. Instead, gravitational collapse of filaments and/or merging 
of filamentary structures can lead to building up dense high-mass star-forming clumps. However, the overall density structure of the cloud has a 
^H , large influence on this process, and requires a better understanding. 

""> ■ Methods. The Rosette molecular cloud, irradiated by the NGC 2244 cluster, is a template region for triggered star-formation, and we investigated 
^^ ' its spatial and density structure by applying a curvelet analysis, a filament-tracing algorithm (DisPerSE), and probability density functions (PDFs) 
on Herschel* column density maps, obtained within the HOBYS key program. 

Results. The analysis reveals not only the filamentary structure of the cloud but also that all known infrared clusters except one lie at junctions of 
filaments, as predicted by turbulence simulations. The PDFs of sub-regions in the cloud show systematic differences. The two UV-exposed regions 
/-s have a double-peaked PDF we interprete as caused by shock compression, while the PDFs of the center and other cloud parts are more complex, 
j.^ ' partly with a power-law tail. A deviation of the log-normal PDF form occurs at A v «9'" for the center, and around 4'" for the other regions. Only 
+7? , the part of the cloud farthest from the Rosette nebula shows a log-normal PDF. 

jrt ■ Conclusions. The deviations of the PDF from the log-normal shape typically associated with low- and high-mass star-forming regions at A v ~3^1 m 

1 ' ' and 8-10"', respectively, are found here within the very same cloud. This shows that there is no fundamental difference in the density structure of 

low- and high-mass star-forming regions. We conclude that star-formation in Rosette - and probably in high-mass star-forming clouds in general - 

t-H . is not globally triggered by the impact of UV-radiation. Moreover, star formation takes place in filaments that arose from the primordial turbulent 

^ . structure built up during the formation of the cloud. Clusters form at filament mergers, but star formation can be locally induced in the direct 

s ! ' interaction zone between an expanding H n-region and the molecular cloud. 

,_j_ ■ Key words, interstellar medium: clouds - individual objects: Rosette 

fj , 1. Introduction spatially highly inhomogeneous molecular cloud has little or no 

f~0 impact on the global evolution of the cloud due to the presence of 

^ . The role of feedback from massive stars has been long discussed s accretion flows. Generally, in a dynamic view of molec- 

. . ; as a dynamical trigger for subsequent star formation (SF). One mar doud . and SF . interactions (e Ballesteros-Paredes et al. 

m > , has to distinguish, however, between the driving scales on which g^ starg form in itationall co ii apsing den se filaments 

^ . the processes work -and their re alive importance Large-scale ^ formed in turbulen ce-generated, interacting shocks. Dale & 

>£j ; turbulence, even on a Galactic scale (Mac Low & Klessen[2004|), BonneU ^j^ d that the most massive clusters should 

d . can be excited by expanding supernova shells, colliding flows bg found at fce junctions offUaments . This confirms the finding 

forming cloud^or globally by any form of accretion (Klessen of Schneider et al (2010a) where a mass flow th h subfila . 

& Hennebelle IMS- Intermediate- to small-scale processes in- mentSi babl med b magnet ic fields, was observed for 

elude radiation, stellar winds, and outflows. Observationally ex- me DR2 1 filament The diction that (massive) c i uste rs form 

panding ionization fronts were proposed to account for frag- ferentiall in the d tential wells where filaments 

mented shells around bubble-like H n-regions (e.g. Zavagno et can be verified observationall and it is within the sc of this 

al I2II3 Shock compression of a simple layer is the imtialidea to do so We chose the Rosette Molecular cloud (RMC) 

of the collect and collapse scenario (Elmegreen & LadaQ2Z3. as a proposed template for 'triggered' SF. A discussion of the 

Recent models (Walch et al. IMB consider now a n initially in- advanta and disadvantages f this premise can be found in 

homogeneous cloud. In particular Dale & BonneU dMU I2HH Schneider et al. (2010b). To characterize the cloud structure, a 

showed that photoiomzin g radiation from an OB-cluster on a curvelet decomposition of column density maps obtained with 

* TT ... cc . , . ... . . . Herschel data is performed, and a filament-tracing algorithm is 

Herschel is an ESA space observatory with science instruments pro- ,• , , • ™, » 

vided by European-led Principal Investigator consortia and with impor- a PP lled t0 . the ma P ( see Arzoumanian et al. 201 1 and Appendix 
tant participation from NASA. B for details). 
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Fig. 1. Left: Curvelet map of the column density, expressed in logarithmic magnitudes of visual extinction (Ay= 0.05 to 7). The 
absolute values, however, depend on the decomposition threshold (»:20% of the intensity are in the curvelets) and are not relevant 
for this study. The filamentary structure becomes apparent and is additionally traced by DisPerSE (white lines). See Appendix B 
for more details and a close-up of the central region. Known infrared-clusters are indicated as blue stars ('PI' for Phelps & Lada 
119971 A-F for Poulton et al. 2008; and Refi for Roman-Zuniga et al. 2008). White stars in the upper right corner indicate the O-stars 
from NGC 2244. Gray triangles are the massive dense cores (MDCs) identified in Motte et al. (2010). Right: A map showing the 
'confidence' (coded from to 1 in color scale) of finding cluster formation sites. This map combines the information of the column 
density map and the locations of filament junctions (see text for details). Known clusters and MDCs are labeled as in the left panel. 
The dashed line indicates the approximate border beyond which (toward the cluster) SF due to the direct compression scenario may 
be at work. 



The Rosette nebula at a distance of 1.6 kpc is illuminated by 
the central OB cluster NGC 2244 located inside a cavity. The 
expanding Hn-region is interacting with a GMC (mass of a 
few 10 5 M Q , Williams et al. 11995b . Photon dominated regions 
(PDRs) are detected not only along the interface but also deep 
inside the molecular cloud (Schneider et al. 11998b . reflecting 
the deep penetration of UV-radiation owing to its inhomoge- 
nous structure. We investigated the density structure of the RMC 
by deriving probability density functions (PDFs) of the column 
densityQ of the whole cloud and subfields in the Rosette. To do 
this, we used observations from Herschel within the HOBYS 
{Herschel imaging survey of OB Young Stellar objects, Motte, 
Zavagno, Bontemps et al. 2010) guaranteed time key program. 
Column density maps were determined from a modified black 
body fit to the wavelengths of PACS and SPIRE (see Appendix 
A) and have an angular resolution of «37". Generally, a PDF 
characterizes the fraction of gas that has a column density N in 
the range [N, N+AN] (e.g., Federrath et al. 120081 



1 expressed in visual extinction A v withN(H 2 )/Av=0.94xl0 21 cm 
mag -1 (Bohlin et al. [T978l 



2. Results and analysis 



2. 1. The filamentary structure of Rosette 

Figure Q] (left) shows a curvelet image from a multi-resolution 
morphological decomposition (Starck et al. I2004I I of the col- 
umn density map of the Rosette (see Fig. 2) obtained from 
Herschel data. Overlaid on the image are the filaments traced 
by the DisPerSE algorithm (Sousbie et al. I201H . This method 
reveals the highly filamentary structure of the RMC, which is 
not fully apparent from visual inspection of the column den- 
sity map (or molecular line maps). Existing infrared (IR) clus- 
ters and the most massive dense cores (MDCs) identified in the 
same data set (Motte et al. 2010; Hennemann et al. ;201O]) are 
overlaid on the image. The latter are potential sites of future 
massive star formation. It is apparent that all sources lie in the 
proximity of junctions of the most prominent filaments in high 
column density regions. Some of the MDCs and the cluster P12, 
however, are also located along filaments (see, e.g., close-up of 
the center region shown in Fig. 5). To better qualitatively char- 
acterize the filament-junction/cluster correlation, we created a 
'confidence map' that combines the two requirements for cluster 
formation: high column density and filament junction. First, we 
produced a mask containing only the junction points and then de- 
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fined a Gaussian distribution function with a FWHIVQ of 3'= 1.5 
pc around each point. This procedure yields a 'junction-area' 
that is more realistic than only a single point. We then multi- 
plied this map with the column density map, and normalized the 
resulting map to obtain values between and 1. Inspection of 
Fig. Q] (right) shows that the IR-clusters and many MDCs are in- 
deed found in regions of high column density where filaments 
merge (all within the maximum cluster size of 7 pc). The only 
exception is the H ii-region/molecular cloud interface (northwest 
of the dashed line in Fig. [TJ. Neither cluster P12 nor the ma- 
jority of MDCs lie in regions of high confidence. Here, shock 
compression of the expanding ionization front with subsequent 
fragmentation into dense cores may be at work, as suggested by 
the PDFs of this region (see next Sec. l2.2l i. 
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Fig. 3. PDFs for the six subregions in the Rosette, labeled from 
1 to 6. The green line indicates the fitted PDF. The Ay -value of 
the deviation of the PDF from the log-normal form is shown as 
a dotted vertical line and the width of the PDF {cr n ) is given in 
each panel. 



Fig. 2. Column density map of the Rosette molecular cloud, ex- 
pressed in visual extinction (see color bar). The six regions that 
were selected to determine individual PDFs are indicated by 
polygons. 



2.2. PDFs from the Herschel column density map of Rosette 

Figure [2] shows the column density map including our divi- 
sion into several subregions, reflecting different morphological 
and physical properties. (1) and (2): interface of H n-region and 
molecular cloud with many UV-exposed pillars, (3): dense, main 
SF region of the RMC, including several IR-clusters (P14/5/6), 
(4) and (5): more diffuse and colder gas with dense, SF clumps 
(P17 in region 4 and P13 in region 5), (6): quiescent and cold 
remote part of the cloud without SF activity. 

The PDF obtained for the whole molecular cloud is dis- 
played in Appendix C (Fig. 6). A log-normal PDF with a width 
of cr n =0.63 is fitted up to Ay -3"' (noise level «l m ). The ex- 
cess at high column densities (Ay =3-20'") follows a power law. 
The slope corresponds to a volumetric, radial density profile 
of n oc r~° 65 (Federrath et al. 1201 U . This does not assume 
that the cloud is a single sphere, moreover it consists of many 
clumps/cores, each with a radial profile leading to a common ex- 



2 The radii of the clusters vary between 1.3 and 3.6 pc (Roman- 
Zuniga et al. 2008), as a conservative value we assumed 1.5 pc. 
However, the resulting map does not change much in a range of ~1 
to ~5 pc. 



ponent a. The exponent is smaller than what is typically found 
for dense cores (-1.5 to -2), suggesting that on large (cloud)- 
scales, turbulence is the dominating process compared to grav- 
ity. 

Figure |3]displays the PDFs for the six subregions, all show- 
ing a log-normal distribution at low Ay and a more complex 
curvature or several peaks at higher Ay. Particularly interesting 
are the PDFs of the UV-exposed regions 1 and 2 (the OB -cluster 
is only 10-15 pc away in projected distance) because there are 
two 'peaks' (at Ay«;3 m and 6 m ). A similar double-peak PDF was 
observed (Schneider et al., priv. comm.) using a column density 
map of RCW120 (Zavagno et al. 2010), a very good example of a 
simple, 'bubble' -like geometry. We interpret the second peak as 
arising from a higher column density component due to gas com- 
pression by the expanding ionization front. This interpretation 
is supported by recent numerical models including turbulence 
and radiation (Tremblin et al., in prep.) that exhibit double-peak 
PDFs depending on the turbulent state of the cloud and the initial 
curvature of the cloud surface. We therefore do not expect to see 
this feature in all UV-illuminated environments. Region 3 is the 
central SF region with the highest column density (up to 70'"). 
The PDF here is broad (cr, ; =0.41), shows two 'breaks' (around 
9'" and 20'"), and has a peak value at A v *5 m (all other PDFs 
peak at »2'"). The PDFs of regions 4 and 5 cover a lower density 
regime and deviate from the log-normal form around Ay *4 m . 
The most quiescent region 6 is the only one displaying a well- 
defined log-normal PDF with the smallest width (0-^=0.19). 
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3. Discussion 

3.1. The density structure of the RMC 

In Csengeri et al. (in prep.), we found that low-mass SF regions 
have a break in the log-normal shaped PDF at Ay ~A m (and pos- 
sibly around 8'"), and high-mass SF regions around Ay ~8™ and 
higher. We now show that these two breaks in Ay are found 
in the PDFs of subregions within one single cloud, the Rosette 
GMC, suggesting that there is no fundamental difference in the 
density structure of low-mass and high-mass SF regions up to 
an Ay of typically 7-10"'. Above this value, the tail(s) in the 
PDFs are associated with SF activity. Namely, if gravity starts 
to play a role during molecular cloud formation, an increasing 
fraction of gas will exceed a certain threshold of column den- 
sity/extinction and form stars. Other studies (e.g., Lada et al. 
|20T01 Heiderman et al. l20T0l Andre et al., 2010) have arrived at 
the same value of around 8 m . In addition, only numerical mod- 
els that include gravity are able to reproduce this threshold and 
the corresponding break in the PDF. For example, Ballesteros- 
Paredes et al. (201 1) showed that PDFs vary during cloud evolu- 
tion. Purely log-normal shapes were found in an initially turbu- 
lent, non-gravitating cloud, while one or more log-normal PDFs 
at low column-densities and a power-law tail for higher values 
were found for later stages of cloud evolution. In addition, the 
slope of the power-law tails varies with time. All PDFs in the 
Rosette have the shape of the later evolutionary states, indicat- 
ing that gravity remains an important process in cloud evolution, 
including star-formation. The impact of UV radiation, however, 
needs to be investigated using hydrodynamic simulations includ- 
ing gravity and radiation. 



3.2. Where do stars form ? 

Our Rosette observations indicate that clusters preferentially 
form at filament junctions but cores (and subsequently stars) can 
also form along filaments (see also other studies of high-mass 
SF regions, e.g., the Nessie-nebula (Jackson et al. 1201 Oi l. W48 
(Nguyen Luong et al. 1201 11 1 or Cygnus (Motte et al. 120071 1). 
In Herschel observations of low-mass SF regions (Andre et al. 
120101 Arzoumanian et al. 1201 U dense cores were found to be 
located along super-critical filaments^ without higher concen- 
trations at the junctions. This points toward a scenario in which 
it is the large potential well of merging filaments (subcritical or 
supercritical) that enables the formation of massive stars within 
a cluster. Conform with numerical simulations (e.g. Klessen & 
Hennebelle 120 1 01 Vazquez-Semadeni et al. 2009) it is therefore 
ultimately the total mass and the high accretion rate that is the 
decisive factor for the formation of a cluster. Observationally, 
apart from the present study, there are more examples of how 
massive structures can be built up from parsec-scale filaments 
(DR21, Schneider et al. I2010al and Hennemann et al., in prep.; 
Vela C, Hill et al. 2011). The supply of material may con- 
tinue on much smaller (subparsec) scales to build-up a cluster 
with massive stars, as proposed by, e.g., Csengeri et al. (1201 1\ . 
Alternatively, massive stars may form in a turbulence-regulated, 
quasi-static scenario (McKee & Tan 120031 ). However, the dis- 
cussion of these small, subparsec scale processes is beyond the 
scope of this paper. 



3.3. Does triggered star-formation exist ? 

Dale & Bonnell (2011, 2012) simulated the gravitational col- 
lapse of a turbulent (giant) molecular cloud exposed to pho- 
toionizing radiation. Their main conclusion was that although 
the gas and dust is heated, photoionization has no strong global, 
i.e. tens of parsec scale, effect on the dense gas because the mas- 
sive clumps out of which OB clusters form are continuously ac- 
creting and most of the photoionizing flux is absorbed by these 
accretion flows. Observationally, no clear picture emerges . In the 
RMC, a very low UV-field in the remote part of the cloud was 
determined (^2-8 Habing at a distance of 30 pc, Schneider et 
al. 1998), supporting the idea of low UV-impact. On the other 
hand, a temperature gradient, and a tentative age-gradient of 
sources, was derived (Schneider et al. 2010bJJ In addition, lo- 
cal triggering of SF, as seen in the RMC (and RCW36, Minier 
et al., in prep.) is possible (Sec. 2.1). The observed double- 
peak PDFs (Sec. 2.2) support this scenario. We therefore con- 
clude that though the morphological and PDF analysis of Rosette 
shows no evidence that SF is globally triggered (across the whole 
extent of the cloud), but is moreover governed by gravity, the im- 
pact of UV-radition is not settled and needs more investigation. 
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Appendix A: Observations and column density map 

The Rosette was observed during the Science Demonstration 
Phase (SDP) within the HOBYS key program. This program is 
dedicated to the earliest stages of high-mass star formation and 
images all molecular complexes that form high-mass stars at less 
than 3 kpc with SPIRE (Griffin et al. l20T0t and PAC S (Pog litch 
et al. 120 lOt using the Herschel satellite (Pilbratt et al. 120101 . The 
SPIRE and PACS dat<0 from 70 fim to 500 p.m were obtained 
on October 20, 2009 in parallel mode with a scanning speed 
of 20"/sec. Two orthogonal coverages of size l°45'xl°25' were 
performed, mapping the (largest) southeast part of the Rosette 
Molecular Cloud. The SPIRE data were reduced with HIPE ver- 
sion 7.1956, using a modified version of the pipeline scripts, i.e., 
observations that were taken during the turnaround at the map 
borders were included, the most recent calibration tree was used, 
and the destriper-module with a polynomial baseline of 0th or- 
der was applied. The two orthogonally scanned maps were then 
combined using the 'naive-mapper' (i.e., a simple averaging al- 
gorithm). The destriper module significantly improved the re- 
sulting maps, compared to the first data reduction just after the 
SDP (Schneider et al. 2010b). The angular resolutions at 160, 
250, 350, and 500 jum, are -12", -18", -25", and -37", re- 
spectively. 

The column density was determined from a pixel-to-pixel 
modified black body fit to four wavelengths of PACS and SPIRE 
(160, 250, 350, 500 //m, all maps were smoothed to the beam- 
size of the 500 //m map, i.e., -37"). For the region covered 
by PACS and SPIRE simultaneously, we fixed the specific dust 
opacity per unit mass (dust+gas) approximated by the power law 
Kv = 0.1 (v/lOOOGHzf cm 2 /g and/3=2, and left the dust tem- 
perature and column density as free parameters (see Hill et al. 
201 1; Arzoumanian et al. 201 1 for details). For the region only 
covered with SPIRBj we used 17.3 K, the median value of the 
SED-derived temperature across the main map where PACS and 
SPIRE overlap giving us four-band coverage at 160, 250, 350, 
and 500 //m. We checked the fitted SEDs and found no major dis- 
crepancy in the fits, though the method described above assumes 
a single temperature and optically thin emission, which is not 
always a good approximation if there is a temperature variation 
along the line-of-sight, noise, and if the dust becomes optically 
thick at shorter wavelengths. Shetty et al. (120091 . for example, 
showed that this can even produce an anti-correlation of /? and T 
which was claimed to be observed by other authors. However, it 
is beyound the scope of this paper to go more into detail. While 
the column density map shown in Schneider et al. (2010b) was 
calibrated using extinction maps, we now recovered the Herschel 
zero-flux levels of the Rosette field with Planck data (Bernard 
et al., priv. communication). The final column density map is 
shown in Fig. 4. 



the curvelets (a difference of around +10%, however, does not 
change the overall picture). This reveals the filamentary struc- 
ture without completely suppressing the more compact sources. 
The DisPerSE algorithm detects filaments starting from a given 
threshold (defined as difference between saddle points and peak 
values) on the curvelet image. However, the column density map 
is a 2D-projection of the volume density while DisPerSE works 
topologically, connecting all emission features such that projec- 
tion effects may create links between filaments that are not phys- 
ically related. To overcome this caveat, filament tracing using 
molecular line data cubes can be a solution, as first tests on the 
DR21 filament have shown (Schneider et al., in prep.). 

Figure|5]shows as an example a close-up of the crowded cen- 
ter region of Rosette where different thresholds were applied. A 
threshold is defined by an intensity contrast between pixels, the 
lowest one starting at 0.5xl0 21 crrr 2 and accordingly finding 
many filaments, up to to a fairly conservative value of 2.9xl0 21 
crrr 2 , leaving only the most prominent features. For this paper, 
a threshold of l.OxlO 21 cm" 2 was choosen in order not to detect 
too faint filaments. It is only the most prominent ones that are 
able to provide enough material that is accreted onto the clus- 
ter center to build up high enough masses. Again, changing the 
threshold does not alter these prominent filaments, they always 
remain detected. 



Appendix C: Probability density function of the 
Rosette cloud 

The probability density function of the whole cloud obtained 
from the column density map (Fig. 2 or 4) is shown in Fig. 6 in 
linear and logarithmic scaling. It displays a log-normal form for 
lower column densities and a clearly defined power-law tail for 
higher column densities that was fitted with a power-law. This is 
not generally the case, while the PDF of the high-mass SF region 
NGC6334 also shows a clear power-law tail (Russeil et al., in 
prep.), the PDFs of other high-mass SF clouds are more complex 
and show several breaks in the PDFs (Hill et al. 201 1, Csengeri 
et al., in prep.). The reason for that can partly be line-of-sight 
effects and limited angular resolution. A more detailed disussion 
of PDFs of intermediate- and high-mass SF regions and compar- 
ison to models is presented in Csengeri et al. (in prep.). 



Appendix B: Curvelet/wavelet decomposition and 
filament tracing 

There are two parameters determining the detection of filaments: 
First, the seperation into curvelets and wavelets, and second 
the threshold of DisPerSE to detect filaments. From our expe- 
rience on the curvelet/wavelet analysis from Herschel column 
density maps so far (Andre et al. 2010; Arzoumanian et al. 201 1; 
Hill et al. 2011), and from several tests with Rosette, we ar- 
rived at a good compromise of 20% of the intensity being in 



3 Data are public available in the Herschel Science Archive. 
6 The instruments are offset by 11' in the focal plan. 
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Fig. 4. Column density map of the Rosette Molecular Cloud, obtained from the 160, 250, 350, and 500 jjm maps from Herschel. 
Known infrared-clusters are indicated in the plot as blue stars ('PI' refers to Phelps & Lada (11997) , A, B, etc. the clusters from 
Poulton et al. ( 2008 ), and labeling "Refl" those from Roman-Zuniga et al. (2008)). White stars in the upper right corner indicate the 
O-stars from NGC2244. 
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Fig. 5. Close-up of the curvelet image (Fig. 1) of the center region of the Rosette molecular cloud with the filamentary structure 
traced by DisPerSE, indicated with white lines. Different thresholds (called 'persistence', i.e. intensity contrast level, for filament 
detection were used (2.9, 1.0, 0.5xl0 21 cm~ 2 from left to right). The curvelet image has a maximum value for the column density 
of 2xl0 22 cirT 2 with a sigma of 0.5xl0 21 and the original column density map a maximum of 9.4xl0 22 ctrT 2 with a sigma of 
1.4xl0 21 . 
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Fig. 6. Probability density function of the whole cloud obtained 
from the column density map in linear (top) and logarithmic 
(bottom) scaling. The upper labeling indicates the visual extinc- 
tion. The red dashed line shows a power-law fit (the high-density 
range beyond Ay =20'" is not well resolved, we therefore do not 
attempt to fit a second power law.) 



